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ABSTRACT 

The temperature distribution of the emitting plasma is a crucial constraint when studying 
the heating of solar flare footpoints. However, determining this for impulsive phase footpoints 
has been difficult in the past due to insufficient spatial resolution to resolve the footpoints from 
the loop structures, and a lack of spectral and temporal coverage. We use the capabilities of 
Hinode/^IS to obtain the first emission measure distributions (EMDs) from impulsive phase 
footpoints in six flares. Observations with good spectral coverage were analysed using a regular- 
ized inversion method to recover the EMDs. We find that the EMDs all share a peak temperature 
of around 8 MK, with lines formed around this temperature having emission measures peaking 
between 10^* and 10^^ cm~^, indicating a substantial presence of plasma at very high temper- 
atures within the footpoints. An EMD gradient of EM(T) ~ T is found in all events. Previous 
theoretical work on emission measure gradients shows this to be consistent with a scenario in 
which the deposited flare energy directly heats only the top layer of the flare chromosphere, while 
deeper layers are heated by conduction. 

Subject headings: Sun: activity - Sun; chromosphere - Sun: flares - Sun: transition region - Sun: UV 
radiation - Sun: X-rays, gamma rays 



1. Introduction 

A clear observational description of the plasma 
properties of the lower atmosphere footpoints of 
solar flares provides a critical constraint on the dis- 
tribution of the flare excess energy in this region, 
and hence the profile of flare energy deposition 
and its possible modes of transport. The emission 
measure distribution (EMD) is a way to describe 
the amount of emitting plasma as a function of 
its temperature, and in this paper we present the 
first determination of impulsive phase flare foot- 
point EMDs made using the Extreme Ultraviolet 
Imaging Spectrometer (EIS, ICulhane et al. 2007 ) 
onboard the Hinode satellite. The temperature 
coverage of EIS makes it extremely well suited to 
studying the properties of footpoints during flares, 
where impulsive stage temperatures can be very 



high - on the order of 8-10 MK, but present in the 
lower atmosphere during the extreme conditions 
of a flare. 

There have been several previous studies of the 
distribution of emission measures of solar flares. 
For example, EUV data from the Skylab NRL 
slitless spectroheliograph enabled spatial as well 
as spectral information to be deduced, though the 
"overlappograms" produced convolved spatial and 
spectral information leading to source confusion. 
Several flare EMDs from the rise and decay phases 
have bee n published from this instrument. The 
EMD of iDere et"all (|l977l) . from just before the 
peak of the 0.5-3 A rise phase of flare SOL1973- 
09-05T18:3^, which would correspond roughly to 



^ we use the so lar observation target id entification convention 
described bv lLeibacher et all ||201(]| 1 
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the end of the impulsive phase, had a steep slope 
up to a peak at logT r-^ 6.9. There was evidence 
for emission at Fe Xiv to Fe xvi at concentrated 
footpoints or near the ends of loops, as well as in 
loops themselves. 

The distribution of emission measuretn found 



in the impulsive phase of a flare by IWidinj ( 1982 ) 
in SOL1973-12-22T00:24 had slope / n^irfdS ~ 
2^0.8 YogT = 5.4 up to logT = 6.2, with 

the suggestion of a decrease above logT ~ 6.9. 
Here, rie is the electron number densit y and S the 
distan ce along the line of sight. iWiding fc Hiei 
(jl984l ) presented emission measure distributions 
from two compact flare sources just after the 
impulsive peak of SOL1974-01-21T23:24, one of 
which had a shallow slope (scaling as T°-^) up 
to a maximum at logT 6, and the other 
had a slope of 3 up to logT 6.9. iWidin3 
(jl982( ) emphasizes that these slopes are differ- 
ent from the 3/2 slopes which are often assumed 
to be characteristic of flares. Other observations 
have been obtained during decay phases of flares 
(jPere fc Cooklll979t IWiding fc Spiced Il98r< ). and 
have also tended to attract mor e theoretical atten- 
tion (jMachado fc Emslie|[l979l ). 

High densities in impulsive-phase sources were 
also found using the Sky lab data, inc l uding 



in high temperatur e lines (|Doschek et al.l 11977 



Feldman et al.l 119771). Related observations from 



the Yohkoh satellite (jActon et al.lll992i) have also 
shown high temperatures and densities during the 
impulsive phase, and more spec ifically soft X-ray 



emiss i on from the f o otpoi nt s (iMcTiernan et al 



1993 : iHudson et al.l 119941 : iMrozek fc Tomczak 
20041 ). High footpoint temperature and density 
observations have recently been amply confirmed 
by EIS (e.g.lMilliganlbOllUDel Zanna et al.l l201ll: 



Graham et al.ll201ll : IWatanabe et al.ll2010[ ). 

In this paper, we are able to isolate the foot- 
points in six flares, and determine their impulsive 
phase EMDs using EIS. The ~ 4" spatial reso- 
lution of EIS is sufficient to make a clear iden- 
tification of footpoints and loops, based on their 



distinction should be made here between the EMD and 
differential emission measure (DEM). The DEM f(T) is 
the quantity of units cm~^K~^ derived from an inversion 
of the data, e.g. via the method described in Section 4. 
Integrating the DEM over a fixed logarithmic temperature 
interval gives the emission measure as a function of tem- 
perature EM{T) in the more practical units of cm""" 



morphology at different wavelengths. We iden- 
tify six flares having consistent spectral cover- 
age and raster observations during their impul- 
sive phases. The EMDs recovered are remarkably 
consistent with one another, peaking at a tem- 
perature of logT ~ 6.9, at a maximum value of 
1028 _ io29cm-5 and with EMDs below the peak 
temperature characterized by EM{T) ~ T. 

Section [2] details the EIS data reduction, while 
an overview of the flare observations and their se- 
lection is found in Section [S] Section 4 describes 
the inputs required in calculating the DEM and 
the regularized inversion method used to deter- 
mine the DEM. The footpoint EMDs are shown 
in Section 5 and discussed in depth in Sections 6 
and 7. 

2. EIS Observations 

To obtain EMDs of flare footpoint plasma we 
require observations of flares during the impul- 
sive (rise) phase at multiple temperatures, ide- 
ally in emission lines given their narrow sensitivity 
to temperature (see Section H]). Its spatial, tem- 
poral, and spectral capabilities make Hinode/FJS 
well suited for this task. The spectrometer slit 
can be scanned across the area of an active re- 
gion multiple times during the flare, building up 
raster "images" in many wavelengths simultane- 
ously. Good datasets for flare footpoint studies 
are hard to find. Positioning the slit over the small 
footpoints (typically 2-5") early in the flare is not 
always possible, and since March 2008 telemetry 
from the Hinode spacecraft has been limited, re- 
stricting the temperature sampling. 

Flare observations using the EIS raster study 
GAM.ARTB_RHESSI_b_2 fulfilled our require- 
ments and six were selected. These use the 2" 
slit to scan a 40" x 140" area in 3 min 52 s. 
Around 30 lines are present in these rasters but 
we have narrowed the selection down to 15 to best 
cover the temperature range (Figure [T]), and avoid 
density-sensitive or optically thick lines. Line de- 
tails are listed in Table [1] and contribution func- 
tions in Figure [T] All of the data has been cali- 
brated and fitted using the standard eis_prep and 
ElS_AUTO_FiT SolarSoft routines. An exception is 
the analysis of Ca xvii which uses the method 
described below. We also correct for a measured 
17" North-South and 1" East- West offset between 
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the two separate wavelength bands on the instru- 
ment's CCD. 

Table 1: Emission lines selected for EMD analysis 
with rest wavelengths and peak formation temper- 
atures. 



fon 


A (A) 


logwT (K) 


V 


248.460 


5.4 


VI 


184.118 


5.5 


Fe VIII 


185.213 


5.7 


Mg VI 


268.991 


5.7 


Si VII 


275.361 


5.8 


Fe X 


184.537 


6.1 


Fe XI 


188.216 


6.2 


Fe XII 


195.119 


6.2 


Fe XIII 


202.044 


6.3 


Fe XIV 


274.204 


6.3 


Fe XV 


284.163 


6.4 


Fe XVI 


262.976 


6.4 


Ca XVII 


192.853 


6.8 


Fe XXIII 


263.766 


7.2 


Fe XXIV 


192.028 


7.2 



We extract for each emission line the fitted in- 
tegrated line intensity averaged over a 2" x 3" 
area around a footpoint (Ix 3 pixels) centred on 
the pixel brightest in Fe viii. The strongest foot- 
point emission in these events mostly appears over 
one slit position 2" wide. To bin further in the x- 
direction would sample too much of the surround- 
ing area, thus our binning accounts for the spatial 
extent of the footpoint and covers some offset be- 
tween wavelengths. 

A pre-flare background was not subtracted from 
the footpoint emission. This is allows us to make 
consistent comparisons with the EMDs described 
in the introduction, and the theoretical work dis- 
cussed in Sectionini which treats emission from the 
entire emitting column, not only the flare excess. 

A number of the lines selected are blended by 
neighboring transitions, however the true inten- 
sities can be recovered using fitting techniques 
and other observed lines within the raster. The 



CHIANTI v7.0 atomic database ( Dere et al.lll997 



Landi et al.l 120121 ) is used to identify lines con- 
tributing to the measured line profiles. 

The low temperature Fe viii 185.213 line has 
two reported blends of Ni xxiv 185.166 and Ni 



XVI 185.230. Fitting with three Gaussians re- 
veals small 10-20% blue wing contributions of Ni 
XXIV 185.166 around, but not necessarily within, 
the footpoints. In the red wing CHIANTI predicts 
Ni XVI 185.230 to be the stronger contribution, 
however, no obvious third Gaussian could be seen 
from our fitting and so Ni XVI should not have a 
significant effect on the Fe viil intensity. 

Fe XXIV 192.028 normally dominates in flare 
loop conditions and contributions from blends 
of Fe VIII and Fe xi are small. This work con- 
centrates on an earlier phase in the flare where 
the Fe XXiv emission will be much fainter, there- 
fore a significant Fe XI 192.021 contribution must 
be removed. We do this via the method de- 



scribed in iDel Zanna et al.l (j201l[ ). The observed 
Fe XI 201.734 intensity forms a known ratio 
with A192.021, this can then be used to esti- 
mate the A192.021 intensity in the footpoint. We 
find the ratio of 192.021/201.734 in background 
moss/active regions to be 0.43. The Fe viii 
192.004 intensity is predicted by CHIANTI to be 
around 10% of the Fe xi line and is almost negli- 
gible. Accounting for these contributions removes 
most of the active-region emission (seen in Figure 
13]) leaving a mean background level lower than 
the line uncertainties (see Section 14.11 on DEM 
uncertainties) . 



10"' 




5.0 5.5 6.0 6.5 7.0 7.5 
logio T [K] 



Fig. 1. — Contribution functions G{T) for emis- 
sion lines calculated with coronal abundances, 
CHIANTI's own ionization equilibrium file, and 
a constant density of lO^^cm^'^ 



3 



Ca XVII 192.853 forms part of a complex group- 
ing of six O V and two Fe xi lines, therefore 
interpretation of the line is challenging but is 



well documented (IKo et al 



2OIOI : Del Zanna et al. 2011 ) a nd we have adopted 



the approach in iKo et al 



2009t lO'Dwver etal 



(|2009[ ) to estimate 
the Ca XVII intensity. First the two Fc xi 
192.813/192.901 intensities are estimated using 
known ratios to Fe xi 188.216. In most of the 
footpoints the strong density sensitive O v 192.904 
line can be resolved and the other five O V lines 
are estimated from CHIANTI v7.0 using a fixed 
density of lO^^cm"'^. 

3. Flare Observations and Footpoint Se- 
lection 

Six events, varying in GOES class from B1.8 
to Cl.l, were selected from a period between 
May and December 2007 (Figure [2|). Two of 
the six events labeled Event (b) (SOL2007- 
12-14T14:16:30) and Event (f) (SOL2007-05- 
22T23:25:50) have been examined in greater de- 
tail; Event (b) by iMilligan fc DennisI (|2009l ) and 
Milliga — on the subjec t of non-thermal 

line br oadening, and Event (f) bv lDel Zanna et aT 
(|2011[ ) . These events are observed to exhibit foot- 



point EUV and HXR emission, chromospheric 
evaporation, and footpoint electron density en- 
hancements. 

Rasters shown in Figure [3] illustrate the flare 
appearance at 500,000 K (Fe v iii) and at 16 MK 
(Fe xx iv). From the work bv IMilligan fc Dennis 
(|2009l ) and lGraham et all (|201lh . compact bright- 
enings in transition region lines can often be as- 
sociated with RHESSI HXR observations reveal- 
ing the flare energy deposition site. Rasters have 
been selected where compact, co-spatial Fe Viii 
and Fe xxiv emission rises dramatically compared 
to the background. Figure [3^ (Event (a)) shows 
small bright sources appearing in Fe Viii during 
the rise phase. Compact, hot Fe xxiv emission is 
also present at this early stage but becomes more 
significant later as evaporating hot material be- 
gins to fill loop structures. Figure [3b, (Event (b)) 
is sampled slightly later in the impulsive phase 
and shows a hot flare loop forming next to the 
footpoint. 

These events are highly impulsive. GOES 
lightcurves for each event in Figure [2] show that 



the rise phase of most events lasts 2-4 minutes, 
with the longest just under 10 minutes. Typi- 
cally HXR observations are used to verify that 
the EUV emission corresponds to the impulsive 
phase. However, RHESSI HXR data were not 
consistently available, and so we systematically 
use the G OES derivative as a proxy for the HXR 
emission ( Neupert 19681 ): a dashed line on the 
lightcurves marks where this peaks and is used 
as a guide to select EIS rasters during the impul- 
sive phase. In each event at least one EIS raster 
was found spanning part or all of the impulsive 
phase. Rasters have been chosen as early as pos- 
sible in the flare whilst still showing strong EUV 
enhancements, and keeping the GOES derivative 
peak within the chosen raster limits (dotted lines 
Figure [2]) . The flare evolves as the spectrometer 
slit scans right to left over the window. The time 
at which the slit crosses a footpoint is marked by 
a diamond on the lightcurve. 

Given this morphology throughout the EIS 
temperature range, plus supporting wide held of 
view imaging fro m the X-Ray Tel escope (XRT) 
onboard Hinode (iGolub et al.l l2007h and TRACE 
( Handv et al. 19991) (omitted here for space con- 
straints), we are confldent in identifying the flare 
footpoints, as marked with white arrows in Fig- 
ure [H 

4. DEM Technique 

The set of line observations, I\, and corre- 
sponding contribution functions, G\{T,ne), are 
related to the differential emission measure of the 
plasma via Ix = j Gx{T,n^)i{T)dT where the 
DEM is defined as £,{T) = nl{dS/dT). Obtain- 
ing ^(T) from I\ is an "ill-posed" inverse problem, 
with uncertainties in the data resulting in the so- 
lution being non-unique. However, using physical 
constraints on the data can help recover a use- 
ful solution. There are a variety of techniques to 



recover the DEM (c.f. iFludra fc Svlweste: 



ques to 
^ 119861: 



Monsignori Fossi fc Landini Il991 ) with a Markov 
chain Monte Carlo (MCMC) approach being rou- 
tinely adopted for spec tral line data inversion 
(jKashvap fc !Draki[i99i). [n t his paper we use a 
regul arized inversion method (jHannah fc Kontar 
20121 ). which is able to produce solutions similar 
to the MCMC approach but is computationally 
quicker and estimates both horizontal and verti- 
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15 1 5:50 1 5:55 1 6:00 1 6:05 16:10 
Stort Time {14-Dec-07 15:41:40) 



(e) 



06:18 06:20 06:22 06:24 06:26 06:28 06:50 
Stort Time (16-Dec-07 06:16:40) 



(f) 



23:20 23:24 23:28 23:32 
Stort Time {22-Moy-07 23:16:40) 



Fig. 2. — GOES l-SA lightcurves of each event. A pair of dotted vertical lines represent the start and end 
times of the EIS raster and a single dashed line marks the GOES derivative maximum for the event. The 
time at which the spectrometer slit scans each footpoint is marked with a diamond on the light curve. 



cal errors on the DEM solution. Regularization 
adds a "smoothness" constraint to the DEM so- 
lution so that a stable inversion can be recov- 
ered, avoiding a mplification of the uncertainties 
(jTikhonovlll963() . However, this may be a strong 
assumption. For example, the DEM of a loop 
in conductive equilibrium theoretically has a dis- 
continuous, high temperature cut-off. The min- 
imum of the EM loci curves (the curves repre- 
senting the isothermal emission in each line) are 
used as an initial guess DEM solution. A multi- 
thermal DEM solution must be below these EM 
loci curves (since the isothermal solution gives the 
maximum possible emission at that temperature) 
and the regularized solutions achieve this (see Fig- 
ure 2]). The approach iterates until a positive 
DEM is found and also minimizes the chi-squared 
between the measured and regularized line inten- 
sities. For a full explanation of this approach on 
EIS line data compared t o the MCMC method see 
Hannah k Kontarl (|2012h . 



4.1. DEM Inputs and Uncertainties 

Determination of the DEM requires three in- 
puts; line intensities, intensity error estimates, 
and the line contribution functions, here calcu- 
lated using CHIANTI v7.0. These are calcu- 



lated using coronal abundances (jFeldman et al 



19921). CHIANTPs ow n ionization equilibrium 
file (jLandi et al.l 120121 ) and a constant density 
of lO ^^cm"^ suitable for f lare footpoints (see 



e^ 



Watanabe erall I2OIOI: iGraham et allboilt 



Milliganll2011l) . The DEM shape should not be 
strongly affected by density variations, as the 
density-sensitive lines have been removed from 
the analysis. We assume a 20% systematic er- 
ror across all intensity measurements to account 
for the absolute calibration uncertainty between 
lines (P. Young, private communication) and this 
is added to the fitting error. 

It is difficult to determine the corr ect elemental 
abundances for use in flare analysis ( Athav 1994t 
Feldman et al In the standard model, foot- 

point material originates low in the chromosphere 
before being heated to coronal temperatures, yet 
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Fig. 3. — EIS rasters for each event in Fe viii at 500,000K and Fe xxiv at 16 MK showing the morphology 
at low and high temperatures. A white arrow marks the footpoint positions chosen for EMD analysis in each 
raster. Black arrows in panel (b) also highlight the regions used in determining a loop top and active region 
EMD (see Section 6). 



this material rises and mixes with existing coro- 
nal material in loops. The choice of abundances 
used to interpret flare spectra is therefore not 
straightforward. Furthermore, low FIP (first ion- 
ization potential) elements are found to be en- 
hanced in coronal material compared to the pho- 
tosphere while high FIP elements are unchanged. 
Only O V and O vi in our analysis are high FIP. 
The choice of ionization equilibrium is similarly 
uncertain. Wc will investigate the effects of vary- 
ing abundance and ionization equilibrium in Sec- 
tion O 



4.2. Assumptions 

Before interpreting spectra with CHIANTI one 
must be aware of the built-in assumptions. The 
G{T) functions are calculated for an optically 
thin plasma in thermal and ionization equilib- 
rium. Clearly from RHESSI HXR observations 
flare footpoint electron spectra have an inherently 
non-Maxwellian component of the electron distri- 
bution, although this can be small c ompared to the 



total e nergy of the distribution; see lKrucker et al 
(j201l[ ). In addition the footpoint is a location 



of intense heating, so it is possible to assume 
that the plasma is out of thermal and ioniza- 
tion equilibrium. However, the high density of 
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14Dec07 (15:53:10) 16Dec07 (06:22:48) 22lu1ay07 (23:21 :38) 




Fig. 4. — DEMs shown in black for each event in emission measure units cm~^ and the uncertainty hmits 
of the solution by a shaded grey area. The colored curves show measured line intensity divided by the 
contribution function indicating the maximum possible emission, i.e. the EM loci curves. 



the emitting region may provide a high enough 
electron collision rate for the footpoint plasma to 
be close to equilibrium. For a plasma at 10^ K 
and Ue = 10"'^^ cm~^ the electron-electron collision 
timescale is r^e = 1-33 x lO"** s, which is probably 
much shorter than the flare heating timescale, so 
the core electrons of the distribution can rapidly 
reach thermal equilibrium. 

Calcu lation of non-equi librium ionization (NEI) 
states bv lBradshaw ( 2009t ) shows that during heat- 
ing the population of a given ion peaks at a higher 
temperature than that at which it would peak 
in equilibrium. This could lead to systematically 
lower temperatures being inferred. Again the ef- 
fect is more pronounced at lower densities as ion- 
ization and recomb ination processes are driven by 
electron collisions. Bradshaw (l2009l) finds that at 



densities of 10 " cm the peaks can shift by up 
to log T = 0.3 K but by lO^^ ^^-3 ^j^g 

difference is negligible. 

Finally we assume that the plasma is optically 
thin. This may at first sight be a more problem- 



atic assumption given the high densities and foot- 
point emission originating from deeper regions of 
the atmosphere, which may be optically thick. We 
estimate the opacity for Fe viii, Fe x and F e xvi 
from the method in Bloomfield et al.l (2002) and 
MilliganI (|201l[ ). The optical depth at line centre. 
To can be expressed as 



-, -, „ in — 14\i' M TliQj 

To 1.16 X 10 Xhj\ — 

V T Hel 



riH Ne 



Neh (1) 



where A is the line wavelength, the oscil- 
lator strength of the transition, and M the mass 
of the ion. riion/nei and nu/Ne are taken from 
the default CHIANTI ionization equilibrium file 
and coronal abundances. Taking a path length 
of 1" and nn = 5 x 10^ cm~^ we find To = 
0.527,0.06 and 0.04 respectively for Fe viii, Fe x 
and Fe xvi. As expected emission lines at higher 
temperatures are largely unaffected by the plasma 
opacity, yet the cooler Fe viii line may be infiu- 
enced. An optical depth of tq — 0.527 corresponds 
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to a drop in transmission of ~ 40% — compared to 
the line uncertainty of at least 20% this could be 
significant. However, if photospheric abundances 
are used this drops to ^ 10%, demonstrating that 
a careful analysis of flare abundances is required 
in future. To fully understand the effects of non- 
equilibrium and radiative transfer in a flare atmo- 
sphere is a serious undertaking which lies beyond 
the scope of this paper. Our results are the first of 
their kind using the best spectroscopy and atomic 
calculations available, but we must be aware of 
these possible shortcomings. 

5. Flare EMDs 

EMDs derived for the six different events us- 
ing the above methods are shown in Figure |4l 
The figures show the emission measure distribu- 
tion EM[T) in units of cm^^ found by integrating 
the DEM over a fixed logarithmic temperature in- 
terval. This returns an EMD of the same form as 
those discussed in the introduction. In all events 
the EMDs in black lines are bounded by the col- 
ored EM loci curves, confirming that the regular- 
ized solutions are below the expected maximum 
emission. 

A shaded region outlines the extent of the EMD 
uncertainty in temperature and EM space. The 
true solution lies within this boundary. Within 
the plotted temperature range the uncertainties in 
temperature are mostly within an order of magni- 
tude or less in EM. At temperatures above lO*" K 
the solutions have a large uncertainty due to the 
broad G{T) response in Fe xxiii and Fe xxiv. 
Also these emission lines are fainter at footpoints 
so have larger fitting errors. Unsaturated, soft X- 
Ray observations from _ffmorfe/XRT would have 
helped constrain this part of the EMD but were 
unavailable. Extending the temperature range 
beyond the limits of the G{T) functions signifi- 
cantly spreads out the errors at the temperature 
limits. The regularization is unable to find solu- 
tions where G(T) is undefined, hence it is therefore 
not possible to make physical conclusions beyond 
these limits. 

All of the footpoint EMDs share a strikingly 
similar profile: increasing with an almost con- 
stant gradient of EM{T) ^ T to a peak around 
log T = 6.9 then falling off quickly at higher tem- 
peratures. The peak temperature suggests a sig- 



nificant presence of plasma at 8 MK in the flare 
footpoints. Peak emission measures vary between 
lO'^^ and 10^^ cm~^ but do not appear strongly 
related to the GOES class, which is perhaps ex- 
pected for a small sample of lower energy events 
sampled at slightly different times in their evolu- 
tion (see Tabic [2]). 

We measure a gradient of EM (T) ~ T^-'^-!±o.27 
between logT — 5.5 — 6.9 in Event (b), where 
the uncertainty is estimated by the maximum and 
minimum gradients allowed within the EMD error 
region. This footpoint EMD is plotted in Figure [5] 
(green region) against a line of gradient EM{T) ^ 
T resembling our event, and EM{T) ^ T^, a 
commonly observed gradient for a transition re- 
gion/low corona atmosphere. 

To verify that the similarity of the footpoint 
EMDs is not an artefact of the regularization 
method used, we have calculated EMDs from both 
active region (AR) and flare loop top (LT) plasma. 
Again using Event (b), EMDs from AR and LT lo- 
cations are shown in Figure [S] in the orange and 
blue regions respectively, with the footpoint shown 
in green. The AR and LT locations on the raster 
are marked for reference by black arrows on Figure 
[Sh. Their EMDs are noticeably different from the 
footpoint EMD, showing that the regularization 
responds well to different plasma temperature dis- 
tributions. The active region EMD peaks between 
log T = 6.3 — 6.5 with a lower emission measure 
and steep high temperature cut-off. The gradient 
in the AR and LT EMDs are steep (> xi) and re- 
markably similar up to logT = 6.3 where the loop 
top EMD becomes shallower and continues rising 
to over 10 MK — the point at which the EMD 
becomes poorly constrained. 

Emission in the footpoint between tempera- 
tures of logT — 5.2 — 6.2 is greater by up to an 
order of magnitude than in the active region or 
loop top. From the conventional understanding 
of flares it is likely this is chromospheric plasma 
in the process of being heated to flare tempera- 
tures. The break in similarity between the LT 
and AR gradients at logT — 6.3 suggests that 
only plasma above this temperature is being evap- 
or ated into the flare loop, ag reeing with the results 
in lMilligan fc Dennis ( 20091 ) where only footpoint 
plasma above logT = 6.2 exhibits evaporative up- 
flows. A full flow velocity analysis will be the sub- 
ject of future work, but it is reassuring to see that 
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such arguments can be made from the results and 
that the behaviour of the EMDs varies in reason- 
able ways across the event. 




6.0 6.5 
log,o Temperature [K] 



Fig. 5. — A comparison of footpoint, loop, and 
active region EMDs in Event (b). The purpose 
of this is to check the regularization response to 
different temperature distributions. Gradients of 
EM{T) - T and EM{T) - tI are added in 
dashed lines. 




6.0 6.5 7.0 

log,o Temperature [K] 



Fig. 6. — Footpoint EMD of Event (b) using pho- 
tospheric and coronal abundanc es and a variety of 
ionization equilibrium theories ( Landi et al.l 12012 ; 
Mazzotta etaLlll998t iBrvans et al.ll2009l) . EMDs 
are shown by the black lines within a shaded er- 
ror boundary. Curves within the lighter shaded 
area are from photospheric abundances and the 
darker from coronal abundances. For each abun- 
dance file three ionization equilibrium theories are 
tested and plotted in different line styles. 



5.1. Varying the abundance and ionization 
equilibrium 

As discussed in Section 14.11 the appropri- 
ate values of abundance in the chromospheric 
plasma, and of the ionization equilibrium, are not 
known, therefore we repeat the analysis using var- 
ious ionization fil es with photospher i c and coro- 
nal abundance s ( Grevesse fc Sauval ( 1998 ) and 



Feldman et al.l (|l992n respectively). EMD curves 



for the footpoint in Event (b) are shown in Figure 
[6] using photospheric (light shading) and coronal 
(dark shading) abundances. The EMD is around 
a factor of 4-5 larger using photospheric abun- 
dances but varies very little in shape; only at 
very low temperatures where oxygen is dominant 
is there any deviation in gradient. Changing the 
abundances therefore only alters the result signifi- 
cantly by a constant factor in EM, i.e. the gradient 
remains EM{T) - T between logT - 5.5 - 6.9. 

Looking again at Figure [5] changing the ioniza- 
tion equilibrium parameters (solid, dashed, and 
dotted line styles) also has very little effect on the 



EMD. Given the large number of lines and small 
variations these parameters produce on the con- 
tribution functions, the final EMD is relatively in- 
sensitive, especially when considering how much 
larger the intensity uncertainties are. Any differ- 
ences due to different ionization equilibria all lie 
within the error boundaries plotted. 

5.2. Column emission measures 

Estimating the emission measure of emitting 
plasma can also be approach ed through the use of 
density-sensitive line ratios. Milligan (|201ll ) used 
five line ratios available in these rasters (Mg vii. 
Si X, Fe XII, Fe xiii, and Fe xiv) to estimate the 
footpoint electron density at various temperatures 
in Event (b), and from these and the observed in- 
tensities, calculate the respective column depths 
of the emitting material. Here we use these mea- 
surements of column depth for an alternate esti- 
mate of the column emission measure at a range 
of temperatures. 

The intensity of a given emission line, /, inte- 
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Table 2: Event EMD parameters sorted by GOES Class. 



Flare ID 


GOES Class 


Event 


Peak Temp 
logw T (K) 


Peak EM 
xlO^^cm-^ 


SOL2007-12-16T06:22:40 


B1.8 


(e) 


6.9 


4.0 


SOL2007-05-22T23:25:50 


B2.7 


(f) 


6.9 


1.0 


SOL2007-12-14T15:22:00 


B8.8 


(c) 


6.9 


8.0 


SOL2007-12-14T15:54:15 


B8.8 


(d) 


6.9 


2.0 


SOL2007-12-14T01:39:20 


B9.6 


(a) 


6.9 


2.0 


SOL2007-12-14T14:16:30 


Cl.l 


(b) 


6.9 


6.0 



grated over the line of sight column depth, S, can 
be expressed as: 



Mg VII Si X Fc XII 



47r/ ^ 0.83 / G{T,ne)nldS. 



(2) 



By assuming the electron density, where rie is ob- 
tained from independent density diagnostic pairs, 
is constant across each pixel, and calculating the 
line contribution function, G{T,ne), at the mea- 
sure d electron density, the column depth is derived 
(see Milliganl ( 201lh for further details). Since 
the column emission measure EMcoi is defined by 
EM col = J nldS, we have combined the electron 
density and column depth measurements to esti- 
mate EMcoi for each diagnostic line pair at the 
footpoint in Event (b). 

The 5 panels in Figure [7] show maps of col- 
umn emission measure for Mg vii, Si X, Fe xii, 
Fe XIII, and Fe xiv. In each of the five rasters, 
higher column emission measures were found at 
the footpoint locations compared to the surround- 
ing active region, even in the cooler Mg VII and 
Si X lines. The footpoint column emission mea- 
sures returned from Figure [7] are between 10^^ — 
lO^^cm"^, in agreement with the regularized in- 
version method. Uncertainties in both the density 
estimates and regularized EMDs make it difficult 
to comment on deviations between the techniques 
of less than an order of magnitude. However, as 
the observed deviations are not larger than this, 
the column emission measures found using line di- 
agnostics do help to reinforce the emission mea- 
sures obtained via the inversion method. 

6. Discussion 

Different theoretical models of energy inputs 
and losses, during flares and in the quiet Sun, pro- 




1 
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Fig. 7. — The reconstructed column emission mea- 
sures in Event (b), found by combining the measured 
electron densities and column depths for the Mg VII, 
Si X, Fe XII, Fe XIII, and Fe XIV lines. The footpoint 
can be seen at (-90" 560"). 



duce EMDs with different slopes, and the slope 
may provide a diagnostic of the energy balance of 
the emitting plasma. During a flare the energy 
balance equation can be very complex, so it is 
normal to make certain simplifying assumptions 
- for example, that the temperature and ioniza- 
tion fraction of the plasma has reached a steady 
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state, that the plasma is gravitationally strati- 
fied, and that the emitting structure has constant 
cross-section. One of the ear hest treatments by 
Shmeleva fc Svrovatskii ( 19731) investigated a flare 
atmosphere spht into two layers. In the upper 
(high temperature) layer, flare energy deposition 
occurred, and was balanced by conductive losses 
to the lower layer. Conduction was assumed to 
dominate radiative output from this layer. In the 
lower layer, which did not receive any direct flare 
heating (e.g. by electrons), the balance was be- 
tween conductive input from above and radiative 
losses. 

Int erpreting Figure 5 in'Shm eleva fc Sv rovatskii 
(|l973l ). the resuhing EMD gradient in a log EM- 
log T plot using a constant pressure assumption 
was ^ 1.2 over a range of a few xlO^ to lO^K. 
(It is ^ 2.2 using a constant density assumption). 
They argue that constant pressure was a valid 
assumption when dealing with the narrow (com- 
pared to pressure scale height) high-temperature 
transition between the flaring chromosphere and 
corona. 

A similar analysis in Widind (1982), looking 
only at the layer where conduction and radiation 
balance, demonstrated a log EM — log T slope of 
1, as we observed . The difference betwe e n this 
result and that of Shmeleva k. Svrovatskii ( 19731 ) 
may lie in the form of the radiative loss function 
adopted. 

The agreement betw e en th e slopes found b y 
Shmeleva fc Svrovatskiil (Il973l) . IWidind (Il982l) . 
and those derived from our observations, is in- 
triguing. Perhaps it suggests that in our events 
energy injected at the footpoints is localized to the 
very top of the flare chromosphere, in a region at a 
temperature of logT > 6.9, with the temperature 
structure beneath determined primarily by con- 
duction and radiation. In other words, any direct 
flare energy input in the region log T ~ 5.5 — 6.9 is 
negligible in magnitu de comp a red to other energy 



loss or gain terms. Brosiud ( 2012 ) has inferred 



similar behaviour in a C6.6 class flare. The early 
appearance of coronal Fe xix emission and late rise 
of transition region lines (O V, Si xii and He i) 
was interpreted as evidence for transition region 
plasma being heated by thermal conduction from 
directly heated coronal plasma. 

However the agreement between model and ob- 
servation may also be coincidental, as there are 



other assumptions in the analyses discussed. For 
example, it is assumed that conduction is deter- 
mined by classical Spitzer conductivity, but the 
strong temperature gradients implied by the small 
vertical extent of the flare transition region mean 
that non-local or saturated flux effe cts may be 
important (e.g. Battaglia et al.l 20091 ). Mechani- 
cal energy loss and enthalpy flux due to plasma 
flows (e.g. evaporation) are also neglected. These 
loss terms may not be impor tant everywhere in 
the E MD temperature range; iMilligan fc Dennis 
(|2009[ ) show that in Event (b) high speed evap- 
orative upflows a r e only present above ~ 2 MK. 



Underwood et al. ( 19781) show that the effect of 



evaporation tends to be to flatten the slope of the 
EMD. Further work is needed to understand the 
effects of relaxing these assumptions, as well as 
on exploring other descriptions for the flare direct 
energy input and the radiative loss function. 

7. Conclusions 

We have obtained the first emission measure 
distribution of the plasma at a fiare footpoint us- 
ing data from Hinode/FlS in conjunction with a 
regularized inversion method. The spectral imag- 
ing capabilities of EIS allows us to separate the 
footpoint EUV spectra - therefore EMD - from 
loop structures; this ambiguity has been a draw- 
back in many previous studies. The resulting 
footpoint EMDs can be characterized by an emis- 
sion measure gradient of EM{T) ~ T between 
log T ^ 5.5 — 6.9 that falls away at higher temper- 
atures, and peak emission measures on the order of 
IQ'28 _ iq29 (2J-Q-5 ^jjg absolutc valuc of emission 
measure in the EMDs is further confirmed by the 
use of density-sensitive line ratios to estimate foot- 
point column emission measures. In previous the- 
oretical work the EMD gradient is found to be sen- 
sitive to the energy transfer methods in a heated 
atmosphere. Our EMD profiles are in rough agree- 
ment with a flaring mechanism depositing energy 
at the top of the flare chromosphere heated to 
logT > 6.9. Deeper layers then radiate the con- 
ductive flux received from the hot layer above. 

Obtaining reliable measurements of footpoint 
EMDs is not only of theoretical interest, as they 
can be used to better estimate synthetic line in- 
tensities in flares. This is useful for identifying 
blends and lines in other instruments such as the 
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SDO/EVE spectrometer, or in calibrating instru- 
mental responses for future instruments (see the 
forthcoming IRIS mission). 

The next step in understanding the heating of 
footpoint plasma is naturally to use RHESSI HXR 
observations to estimate the energy deposited in 
the chromosphere. How does the EMD depend on 
the HXR spectra of the event, are they as consis- 
tent as the EMD themselves? Combined with the 
temperature response in EUV lines a fuller pic- 
ture would be available to compare various flare 
heating models. 
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